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Abstract 

A  model  has  been  designed  and  constructed  for  the  all-liquid,  sodium  borohydride/hydrogen  peroxide  fuel  cell  under  development  at  the 
University  of  Illinois  at  Urbana-Champaign.  The  electrochemical  behavior,  momentum  balance,  and  mass  balance  effects  within  the  fuel  cell  are 
modeled  using  the  Butler- Volmer  equations,  Darcy’s  law,  and  Fick’s  law,  respectively,  within  a  finite  element  modeling  platform.  The  simulations 
performed  with  the  model  indicate  that  an  optimal  physical  design  of  the  fuel  cell’s  flow  channel  land  area  or  current  collector  exists  when 
considering  the  pressure  differential  between  channels,  and  the  diffusion  layer  permeability  and  conductivity.  If  properties  of  the  diffusion  layer 
are  known,  the  model  is  an  effective  method  of  improving  the  fuel  cell  design  in  order  to  achieve  higher  power  density. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

A  novel  sodium  borohydride/hydrogen  peroxide  (NaBFU/ 
H2O2)  fuel  cell  is  under  development  at  the  University  of 
Illinois  at  Urbana-Champaign.  The  cell  has  distinct  advantages 
for  space  applications  due  to  the  all-liquid  nature  and  high  energy 
density  of  the  reactants.  The  NaBFL/FLC^  fuel  cell  group  at  the 
University  of  Illinois  has  been  working  to  steadily  improve  the 
design  and  capabilities  of  the  cell  [1-3].  Due  to  the  novel  nature 
of  the  fuel  cell,  there  is  a  lack  of  modeling  research  into  its 
characteristics,  motivating  this  study. 

The  development  of  a  model  for  the  NaBFL/FLC^  fuel  cell  is 
intended  with  the  goal  of  understanding  some  of  the  characteris¬ 
tics  that  are  important  in  the  design.  Specifically,  the  geometries 
of  the  fuel  flow  channels  and  the  characteristics  of  the  diffusion 
layer  are  considered.  It  is  hoped  that  the  cell  performance  and 
power  output  can  be  improved  for  later  optimization  by  varying 
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the  geometries  and  physical  characteristics  of  the  cell  and  the 
diffusion  layer.  In  this  study,  flow  rates  (indicative  of  pressure 
differentials),  permeability,  conductivity,  and  basic  flow  channel 
geometry  are  varied  [4] . 

For  reference  throughout  the  paper,  the  balanced  reaction 
equations  for  the  NaBE^/F^CE  fuel  cell  are  shown  in  Eqs.  (1) 
and  (2). 

NaBH4  +  4H20  NaB02  2H20  +  8H+  +  8e“  (1) 

H202  +  2H+  +  2e_  — »■  2H20  (2) 

2.  Description  of  model 

The  model  developed  here  employ  the  equations  and  geome¬ 
tries  similar  to  those  used  in  prior  models  [4-6].  As  in  those 
models,  the  system  is  set  up  with  a  series  of  differential  equa¬ 
tions  to  account  for  the  various  interactions  occurring  in  the 
ME  A;  however,  specific  equations  are  adapted  for  use  with  the 
physics  packages  in  the  finite  element  modeling  software  [6]. 
The  MEA  consists  of  a  diffusion  layer,  Nafion  membrane,  and 
another  diffusion  layer.  Appropriate  catalysts  are  deposited  on 
the  diffusion  layers.  First,  the  physical  layout  of  the  model  used 
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Nomenclature 

Electrical  and  Butler-Volmer  terms 

a 

conductivity 

0 

electric  potential 

i 

current 

k 

exchange  current  density 

n 

moles  of  e- 

p 

charge  transfer  coefficient 

F 

Faraday’s  constant 

T) 

overpotential 

R 

universal  gas  constant 

T 

temperature 

w* 

species  concentration 

w0 

species  input  concentration 

a 

reaction  constant 

f 

temperature  constant 

* 

a 

anode  subscript 

* 

c 

cathode  subscript 

*«n 

anodic  diffusion  layer  subscript 

*d2 

cathodic  diffusion  layer  subscript 

* 

m 

membrane  layer  subscript 

*eq 

equilibrium  superscript 

Mass  transport  terms 

Di 

species  diffusion  coefficient 

R 

reaction  rate 

u 

fluid  velocity 

No 

inward  flux 

n 

normal  vector 

N 

species  input  concentration 

I 

species  concentration 

Momentum  balance  terms 

v' 

mean  fluid  velocity 

k(tc) 

permeability 

viscosity 

A p  (pdiff)  pressure  differential 

is  described,  followed  by  a  description  of  the  equations  mode¬ 
ling  the  physical  phenomenon  occurring  in  each  component  of 
the  MEA. 

An  efficient  method  for  modeling  a  large  system  is  to  look 
for  repetitive  structures  and  then  model  one  of  the  smaller  struc¬ 
tures.  This  allows  for  the  creation  of  simpler  models,  but  care 
must  be  taken  to  ensure  the  boundaries  of  the  structures  are  ade¬ 
quately  considered.  A  drawback  of  modeling  in  this  way  is  that 
the  edge-effects  of  the  larger  system  are  not  modeled.  For  PEM 
fuel  cells  with  serpentine  flow  channels,  the  edge-effects  are  not 
as  critical  as  the  reaction  rates  in  the  repetitive  array  of  flow 
channels,  which  contribute  the  most  to  the  output  power  of  the 
cell.  The  repetitive  cross  section  chosen  is  from  the  midpoint  of 
one  flow  channel  to  the  midpoint  of  the  next  and  extends  across 


Note:  The  *  indicates  a  symbol,  superscript,  or  subscript  that  occurs  multiple 
times,  either  as  different  species  or  on  different  terms. 


the  two  diffusion  layers  and  the  membrane.  The  flow  in  the  flow 
channels  are  considered  to  be  at  constant  concentration.  This 
is  a  simplification,  but  reasonably  accurate  since,  due  to  high 
flow  rates  and  concentrated  reactant,  there  is  little  concentration 
change  between  adjacent  flow  channels.  For  visualization  pur¬ 
poses,  Fig.  1  shows  the  cross  section  of  the  fuel  cell  that  the 
model  considers;  each  physical  area  is  labeled.  In  the  expanded 
repetitive  section,  the  black  outline  indicates  the  portion  that  is 
modeled;  thus,  the  flow  channels  are  considered  boundary  condi¬ 
tions  in  the  model.  Portions  of  the  graphite  plate  are  part  of  the 
model  to  incorporate  current  flow  through  this  region.  Finally, 
the  two  diffusion  layers  and  Nation  membrane  are  incorporated. 

The  model  is  a  2D  approximation  of  the  fuel  cell,  and  is 
created  and  simulated  using  COMSOF  Multiphysics.  Once  the 
model  is  physically  created,  it  is  described  using  built-in  or 
user-defined  differential  equations.  To  describe  the  model,  each 
region  has  coefficients  for  every  differential  equation  (the  coef¬ 
ficients  may  be  zero  if  that  particular  equation  does  not  interact 
in  that  region);  also,  the  boundary  conditions  are  described  for 
each  region  where  a  differential  equation  is  applied.  The  soft¬ 
ware  has  built-in  applications  that  allow  the  regions  to  be  easily 
described.  For  example,  a  simple  DC  current  equation  allows 
for  constant  voltage  or  current  boundary  conditions,  while  also 
indicating  the  resistance  in  the  region.  More  complicated  depen¬ 
dencies  can  be  set  up  by  having  the  resistance  be  a  function  of 
some  other  element,  such  as  temperature,  although  temperature 
was  not  considered  in  this  model. 

The  software  then  generates  a  2D  triangular  mesh  contai¬ 
ning  thousands  of  elements.  Using  finite  element  analysis,  the 
software  converges  to  a  solution  through  successive  iterations. 
Finite  element  modeling  allows  for  models  with  thousands  of 
individual  triangular  sections  and  points  to  be  efficiently  solved 
on  a  desktop  computer.  In  finite  element  analysis,  finer  meshes 
converge  to  more  precise  results.  For  this  simulation,  it  was 
found  that  a  model  with  approximately  10,000-13,000  elements 
was  very  precise  (after  modeling  with  100,000  or  more  elements, 
the  solution  varies  only  slightly,  less  than  1  %).  The  lower  number 
of  elements  allowed  for  a  fast  solution  time  and  many  different 
tests  and  variables  to  be  considered. 

2.7.  Electrical  model 

All  the  current  generated  in  the  fuel  cell  must  pass  through 
the  diffusion  layers  and  graphite  to  escape  the  cell.  Sometimes 
an  afterthought  in  the  design  of  the  fuel  cell,  voltage  drops  due 
to  resistive  losses  can  significantly  affect  the  operating  charac¬ 
teristics  of  the  cell.  Fortunately,  modeling  the  DC  currents  in  the 
fuel  cell  is  straightforward  with  the  present  model.  Use  is  made 
of  the  electrical  conduction  equation,  shown  in  Eq.  (3),  which 
has  previously  been  used  to  model  the  currents  in  fuel  cells  [7] . 

V(-crV0)  =  0  (3) 

The  boundary  layer  at  the  diffusion  layer  to  membrane  interface 
is  also  an  important  consideration.  Normally,  voltages  must  be 
equal  at  the  boundaries  between  the  sections,  but  in  the  case  of 
the  fuel  cell,  current  is  generated  at  the  interface.  The  current 
source  occurs  here  because  the  catalyst  layer  is  treated  as  if  it  is 
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Cross-Section  of  PEM  Fuel  Cell 


Fig.  1.  MEA  cross  section. 


located  on  the  surface  of  the  Nation  membrane.  This  is  another 
simplification  since  the  catalyst  in  the  actual  cell  is  distributed 
over  the  diffusion  layer.  While  some  error  is  introduced,  follo¬ 
wing  normalization  of  the  exchange  current  [4],  any  error  caused 
is  thought  to  be  minimal.  The  Butler- Volmer  equation  is  then 
used  to  model  the  current  generated  at  the  membrane  interface 
[8]. 


i  =  io 


exp 


(^) 


—  exp 


(4) 


Eq.  (4),  the  generic  Butler- Volmer  equation,  is  separated  into 
anode  and  cathode  terms  (the  two  sides  of  the  membrane)  and 
includes  the  normalized  reactant  concentrations.  In  addition,  the 
overpotential  of  the  anode  and  cathode  must  be  considered.  Here 
the  overpotential  is  simply  labeled  77,  but  it  must  also  be  separated 
into  anodic  and  cathodic  terms  in  Eqs.  (5)  and  (6)  [5]. 

u;NaBH4 

4  =  *0,a  XT  nl,  (exp(o;a/?7a)  -  exp(-aa/??a)) 
u;NaBH4o 

?7a  =  0dl  —  0m  —  0a^  (5) 


UAH2O2 

i C  =  *0,c  T  T  ^  (exp(ac/?7c)  -  exp(— ctfc/V/c)) 

1CH2O20 

ric  =  4>d2  -4>m~  0cq  (6) 


is  used  to  eliminate  the  exponentials,  forming  Eq.  (8). 
coshW  =  i(exp<*)-exp(-*))  (7) 


u;NaBH4 

/a  =  2/o,a  -p.  T  T 

u;NaBH4o 


COSh(cya/^a) 


ic 


—  2/q,c 


u;H202 

m;H202o 


cosh(ac/>?c ) 


(8) 


Since  the  current  generation  is  reliant  upon  reactants  being  deli¬ 
vered  to  the  membrane,  convection  and  diffusion  in  the  diffusion 
layer  are  crucial  portions  of  the  model.  At  the  same  time  reac¬ 
tant  depletion  and  product  generation  occurring  at  the  membrane 
interface  depends  on  the  current. 

The  electrical  boundary  conditions  are  very  simple;  the  two 
graphite  ends  of  the  cell  are  considered  to  be  at  constant  vol¬ 
tage.  One  side  is  the  fuel  cell  voltage,  the  other  is  0  V.  In  the 
equations,  it  is  only  the  difference  in  the  two  ends  that  affects 
the  results  so  this  could  represent  any  cell  in  a  fuel  cell  stack. 
Internal  boundary  conditions  simply  follow  Kirchhoff ’s  Circuit 
Laws.  The  top  and  bottom  boundary  conditions  are  electrically 
floating  since  the  repetitive  array  can  be  mirrored  along  these 
lines,  meaning  current  does  not  flow  [6]. 


The  Butler- Volmer  equation  can  be  easily  modified  to  a  more 
useful,  alternative  version  that  appears  to  be  more  robust  when 
solving  the  equation  numerically  while  also  being  mathemati¬ 
cally  identical.  In  this  version,  the  hyperbolic  identity  of  Eq.  (7) 


2.2.  Fluid  model 

There  are  many  different  ways  to  address  momentum  balance 
in  a  fluid  velocity  field,  but  the  diffusion  layer  in  a  fuel  cell 
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is  relatively  unique  due  to  its  porous  nature.  For  this  reason 
it  is  convenient  to  use  Darcy’s  empirical  law  for  laminar  flow 
in  packed  beds,  indicated  in  Eq.  (9).  In  Darcy’s  law,  the  fluid 
velocity  is  proportional  to  the  change  in  pressure  and  inversely 
proportional  to  the  viscosity  [9]. 


The  velocity  profile  and  pressure  differential  are  solved  for  in 
the  model,  but  the  permeability  and  viscosity  are  dependent  on 
the  fluids  and  the  diffusion  layer  used. 

The  mass  balance  in  the  fuel  cell  incorporates  the  fluid 
velocity  from  the  momentum  balance  with  the  movement  of 
the  reactants  and  products.  Modeling  mass  balance  in  a  two- 
phase,  hydrogen/oxygen  (H2/O2)  fuel  cell  is  quite  complicated. 
The  reactants  arrive  as  gasses  but  also  must  diffuse  through 
a  water  layer  at  the  membrane.  However,  for  the  all-liquid 
(i.e.,  single-phase)  NaBH4/H202  fuel  cell,  mass  balance  is 
greatly  simplified.  The  single-phase  system  is  modeled  in  a 
straight  forward  manner  using  standard  convection  and  diffu¬ 
sion  equations.  These  are  even  simpler  if  temperature  gradients 
are  considered  to  be  minimal,  which  is  a  reasonable  assump¬ 
tion  in  the  present  case  given  steady- state  operation  of  the 
fuel  cell  and  the  high  thermal  conductivity  of  the  liquid 
reactants.  The  standard  diffusion  equation  is  Fick’s  Law  in 
Eq.  (10). 

V  •  (- Di Vi)  =  R-u  Wi  (10) 

The  subscript  and  symbol  i  refers  to  the  species,  or  reactant 
and  product  in  the  region.  For  example,  the  anode  has  sodium 
borohydride,  sodium  metaborate,  and  water  for  species. 

Eq.  (10)  must  to  be  considered  for  each  reactant  and  product 
species  used  and  needs  to  be  constrained  to  the  corresponding 
diffusion  layer.  For  the  current  model,  the  reactions  do  not  occur 
in  the  diffusion  layer  itself,  but  instead  only  at  the  membrane 
interface.  Hence,  R  =  0  in  the  diffusion  layers.  The  fluid  velocity 
u  is  determined  using  Darcy’s  law,  requiring  the  coupling  of 
these  differential  equations. 

The  mass  balance  boundary  layers  are  dependent  on  the 
sources  and  the  sinks  in  the  system.  For  reactants,  the  sources 
occur  at  the  inlets  and  outlets  of  the  flow  channels,  which  are 
considered  to  be  at  constant  concentration.  The  sinks  for  the 
reactants  occur  at  the  membrane  interface  and  are  dependent 
on  the  current.  The  sources  for  the  products  occur  at  the  mem¬ 
brane  interface  and  are  also  dependent  on  the  system  current. 
The  flow  channels  are  assumed  to  have  no  products  and  thus 
serve  as  sinks.  The  boundary  layer  equations  for  the  reactants 
have  the  form  shown  in  Eq.  (11). 

-n-N  =  N0  N  =  -DiVi  +  iu  (11) 

In  summary,  the  currents  are  determined  using  the  Butler- 
Volmer  equations,  the  fluid  velocity  with  Darcy’s  Law,  the 
voltages  with  the  DC  media  equations,  and  the  reactant  concen¬ 
trations  with  convection  and  diffusion. 


Table  1 


Fuel  cell  parameters 


Parameter 

Value 

Parameter 

Value 

0  Nation 

15  Sm-1 

Pin 

1.013e— 5Pa 

^Diffusion 

2500  Sm"1 

Pdiff 

500  Pa 

^Graphite 

16670  Sm-1 

£>h2o2 

3.47e— 9  m2  s-1 

K 

1.22e— 11m2 

T>NaBH4 

3e— 9  m2  s-1 

Ma 

1.5  cP 

^NaB02 

1.23e— 9  m2  s-1 

Me 

lcP 

Drag 

3 

3.  Calculations  of  parameters 

Many  of  the  parameters  necessary  to  simulate  the  fuel  cell 
have  been  acquired  through  experimental  means  or  published 
values,  including  the  conductivities,  permeabilities,  diffusion 
coefficients,  and  viscosities.  These  parameters  are  included  in 
Table  1. 

It  is  commonly  accepted  in  literature  that  three  water  mole¬ 
cules  are  dragged  across  the  membrane  with  each  hydrogen 
ion  [8].  We  assume  that  Nafion  drag  coefficient  will  remain 
unchanged  in  the  present  all-liquid  case. 

It  is  also  necessary  to  determine  the  parameters  used  in  the 
Butler- Volmer  equations.  A  Hydrogen  half-cell  was  constructed 
and  used  to  determine  the  exchange  current  density  and  additio¬ 
nal  parameters  such  as  the  Tafel  slope.  The  reversible  potential 
of  each  cell  half  was  determined  using  the  Gibb’s  Free  Energies 
applied  to  the  reactants  and  products  in  each  reaction. 

4.  Simulation  and  results 

Fig.  2  is  the  I-V curve  for  the  model  created.  The  figure  shows 
the  effects  of  the  anode  and  cathode  activation  overpotentials  and 
the  ohmic  losses.  The  cell  achieves  a  power  density  of  about 
0.31  W  cm-2  for  this  region  of  the  membrane.  Some  advanced 
NaBH4/H2C>2  experimental  cells  have  achieved  power  densities 
of  1 .5  W  cm-2,  but  the  normal  version  this  model  treats  typically 
generate  <0.5  W  cm-2  [1-3]. 


I-V/Power  Curve 
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Fig.  3.  Varying  channel  distances  in  serpentine  design. 

Based  on  Fig.  2,  the  model  represents  the  all-liquid  fuel  cell 
fairly  well.  However,  simply  generating  a  model  is  not  the  only 
goal.  Dimensions  and  parameters  within  the  model  are  to  be 
varied  to  better  understand  how  they  affect  the  overall  perfor¬ 
mance  of  the  cell.  First,  the  pressure  differential  between  the  two 
channels  is  varied.  Then,  the  current  collector  land  area  to  chan¬ 
nel  width  ratio  is  varied  versus  the  permeability  and  conductivity 
of  the  diffusion  layer.  The  comparisons  are  done  using  the  I-V 
curves  or  power  density  comparisons,  but  additional  information 
using  other  plots  is  provided  when  useful. 

The  pressure  differential  between  the  channels  is  controlled 
by  two  separate  and  independent  factors:  the  channel  velo¬ 
city  and  the  location  in  the  cell.  An  increased  channel  velocity 
creates  a  higher-pressure  drop  across  the  entire  length  of  the 
flow  channel,  resulting  in  an  increased  pressure  differential 
between  two  neighboring  channels.  The  effects  of  varying  chan¬ 
nel  velocity  (and  hence  a  varying  pressure  differential)  are 
important  to  understand  and  characterize.  In  addition,  the  pres¬ 
sure  differential  is  dependent  on  the  location  within  the  cell. 
The  varying  channel  distance  in  serpentine  designs,  illustra¬ 
ted  schematically  in  Fig.  3,  occurs  because  of  the  serpentine 
design  of  the  flow  channel.  When  the  channel  reverses  direc¬ 
tions  on  itself,  there  will  be  less  pressure  differential  between 
the  two  channels;  the  fluid  has  traveled  less  distance  so  there 
will  be  a  smaller  pressure  drop  between  the  two  points  in  the 
channel. 

The  model  results  shown  in  Fig.  4  indicate  that  an  increa¬ 
sed  pressure  differential  across  the  inlet  channels  raises  the  I-V 
curve,  resulting  in  higher  power.  The  important  characteristic 
of  the  plot  occurs  for  the  pressure  differential  is  at  and  below 
500  Pa;  the  I-V  curve  starts  sloping  downwards  at  higher  cur¬ 
rents.  The  reduced  currents  are  from  mass  transport  effects;  the 
reactants  are  unable  to  permeate  through  the  diffusion  layer  to 
the  membrane  surface  for  reaction.  Intuitively,  this  makes  sense; 
the  higher-pressure  differential  creates  a  higher  fluid  velocity 
between  the  channels,  allowing  more  reactants  to  permeate  the 
diffusion  layer  and  reach  the  membrane. 


Varying  Pressure  Differential  I-V  Curves 


With  the  higher  pressures,  the  fluid  velocity  is  dramatically 
increased  between  the  two  channels  and  reactant  maintains 
higher  concentrations  under  the  current  collector.  The  reactants 
maintain  the  initial  concentration  much  further  under  the  cur¬ 
rent  collector  for  the  higher  velocities.  For  the  remainder  of 
the  simulations,  the  pressure  differential  across  the  channels 
remains  constant  at  500  Pa. 

Next,  the  current  collector  land  area  is  varied.  Basically, 
varying  the  land  area  is  akin  to  varying  the  channel  width,  so 
instead  of  simply  varying  the  land  area,  the  land  area  to  channel 
width  ratio  is  varied.  The  total  width  of  land  area  plus  channel 
is  kept  constant  at  4.4  mm.  By  varying  the  ratio  while  also  using 
multiple  values  for  the  conductivity  and  the  permeability  of  the 
diffusion  layer,  it  is  possible  to  develop  trends  associated  with 
the  land  area.  Fig.  5  shows  the  maximum  power  density  achie¬ 
ved  by  the  cell  versus  the  land  area  to  channel  width  ratio;  as 
indicated,  the  diffusion  layer  permeability  value  is  also  varied, 
in  this  plot. 


Land  Area  Width  vs.  Permeability  (k) 


Fig.  5.  Land  area  width  vs.  permeability  values. 
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There  are  several  competing  effects  illustrated  in 
Fig.  5.  First,  for  the  mid-to-low  permeability  values 
(Kappa <  3.045  x  10“ 12  m2),  high  ratios  greatly  reduce 
the  power  capabilities  as  a  result  of  mass  transport  effects.  As 
discussed  previously,  the  fluid  velocity  under  the  land  area  is 
highly  dependent  on  the  permeability.  With  this  in  mind,  a  large 
land  area  would  logically  create  problems  with  fluid  velocity, 
resulting  in  a  lack  of  reactants  under  the  membrane. 

The  other  effect  seen  in  Fig.  5  is  the  limited  conductivity  of  the 
diffusion  layer.  All  the  current  that  passes  through  the  graphite 
current  collector  must  also  spread  out  to  cover  the  membrane 
surface  to  react  evenly  across  the  entire  surface;  for  narrow  cur¬ 
rent  collectors  (low  ratios),  the  diffusion  layer  provides  enough 
resistance  to  inhibit  this  “spreading  out”  of  the  current  to  cover 
the  membrane  surface.  Fig.  6  illustrates  the  maximum  power 
density  achieved  while  varying  the  land  area  to  channel  width 
ratio  for  several  different  conductivity  values. 

As  the  conductivity  increases,  a  one-to-one  ratio  of  the  current 
collector  land  area  to  flow  channel  width  becomes  optimum. 
For  lower  conductivities,  having  a  wider  land  area  and  narrower 
channel  is  shown  to  be  more  effective.  One  noticeable  trend  is 
the  reduced  power  density  for  all  conductivities  when  the  ratio  is 
less  than  one,  which  illustrates  how  difficult  the  “spreading  out” 
of  the  current  is,  even  at  high  diffusion  layer  conductivities.  All 
the  computations  above  were  conducted  with  a  diffusion  layer 
thickness  of  0.2  mm. 

Different  diffusion  layer  thicknesses  would  probably  result  in 
different  peaks,  but  the  trends  would  still  be  the  same.  The  reason 
for  the  reduction  at  high  land  area  to  channel  width  ratios  most 
likely  results  from  mass  transport  effects;  the  cause  and  effects 
of  mass  transport  have  already  been  discussed  while  varying  the 
pressure  differential  and  the  permeability  values. 

As  the  conductivity  increases,  the  current  achieved  at  the 
center  of  the  flow  channel  increases,  which  intuitively  makes 
sense;  the  higher  conductivity  allows  the  current  to  “spread  out” 
more,  achieving  the  higher  current  density  in  that  region.  The 
narrower  land  area  lowers  the  current  density  in  the  center  of  the 


Land  Area  Ratio  vs.  Conductivity  (a) 


Fig.  6.  Land  area  width  vs.  conductivity. 


flow  channel  for  both  conductivities.  The  narrowest  land  area 
tested  actually  reduced  the  current  density  in  the  center  of  the 
flow  channel  from  13  to  43%,  depending  on  the  conductivity. 
This  indicates  the  difficulty  that  the  current  has  in  reaching  the 
membrane  under  the  flow  channel  due  to  the  increased  resistance 
encountered. 

5.  Discussion 

There  are  several  conclusions  to  be  made  based  on  the  trends 
observed  while  varying  the  pressure  differential  across  the  two 
flow  channels.  First,  operating  at  a  higher  flow  velocity  will 
create  higher-pressure  differentials  throughout  the  fuel  cell, 
increasing  the  reactant  flow  to  the  membrane  and  increasing 
the  average  power  through  the  cell.  However,  high  velocities 
and  the  resulting  pressure  differential  from  inlet  to  outlet  of  the 
reactants  may  not  be  practical.  Second,  the  power  density  in  the 
cell  is  strongly  dependent  on  location  in  the  cell.  As  discussed 
previously,  there  is  a  varying  pressure  differential  across  the  two 
channels  depending  on  the  location  along  the  serpentine  chan¬ 
nel;  this  varying  pressure  differential  creates  a  varying  power 
density  that  is  low  when  the  serpentine  wraps  around  on  itself 
and  high  when  turning  away.  Achieving  a  reduction  in  power 
density  variations  by  redesigning  the  flow  channel  layout  is  not 
a  trivial  task,  and  has  not  been  attempted  though  it  merits  further 
study. 

Finally,  varying  the  current  collector  land  area  to  channel 
width  ratio  resulted  in  observations  as  well.  The  optimum  ratio 
is  highly  dependent  on  the  permeability  and  conductivity  of 
the  diffusion  layer.  The  trends  were  determined  for  the  two 
parameters  independently,  but  in  general,  lower  ratios  reduce 
the  maximum  power  reached  regardless  of  the  permeability  or 
conductivity.  For  low  permeability  values,  the  maximum  power 
density  was  also  reduced  for  large  ratios  due  to  mass  transport 
effects.  For  conductive  membranes,  the  optimum  land  area  to 
channel  width  ratio  is  approximately  one-to-one,  meaning  equal 
widths  for  both.  For  lower  conductivity  membranes,  the  opti¬ 
mum  ratio  increases  slightly.  The  ratio  with  the  highest  power 
depends  on  the  permeability  and  conductivity  of  the  diffusion 
layer,  leading  to  the  conclusion  that  an  optimum  solution  can 
be  determined  through  simulation  if  a  specific  conductivity  and 
permeability  are  provided. 

6.  Conclusion 

A  cross-sectional,  all-liquid  fuel  cell  model  was  created  and 
simulated  using  finite  element  modeling  software.  The  pres¬ 
sure  differential  between  the  flow  channels  and  the  land  area 
to  channel  width  ratio  were  altered  to  determine  the  effects  on 
the  power  density  in  the  model.  Furthermore,  permeability  and 
conductivity  values  were  varied,  with  a  trend  indicating  an  opti¬ 
mum  geometry  exists  to  maximize  the  power  density  for  a  given 
ME  A  permeability  and  conductivity.  Based  on  research  from 
the  studies  performed,  an  optimum  land  area  to  channel  width 
ratio  can  be  determined  for  a  specific  fuel  cell  (assuming  known 
permeability,  conductivity,  and  thickness  of  the  diffusion  layer). 
In  further  studies,  given  geometry  of  the  fuel  cell,  the  optimum 
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characteristics  of  the  diffusion  layer  can  be  determined  through 
simulation,  which  would  aid  the  choice  of  materials. 
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